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Recent rapid advances in the observations of Gamma Ray Bursts (GRBs) lend much 
confidence to the cosmological fireball (FB) model: there is now compelling evidence that the 
radiation is emitted by a relativistic shock where a high speed upstream flow terminates. The 
question concerning what generates such an outflow is central to our search for the ultimate 
trigger mechanism. A key requirement, not well explained by current theories, is that the 
flow must have high entropy-to-baryon ratio. In this Letter we point out that a quantum 
discharge induced by the radiation of the initial FB may be the explanation. The effect is 
likely to be relevant, because the FB energy density inferred from ORB data is large enough 
that the radiation pressure leads to the formation of a surface electric field which is unstable 
to pair creation. Under suitable conditions, such as those of a supernova core, the discharge 
can convert a substantial fraction of the FB energy into surface pairs. This pair plasma is 
not contaminated by FB baryons because it is formed outside the FB. We demonstrate that 
the pairs can then develop relativistic bulk expansion, reaching a maximum speed that meets 
the constraints required to form a GRB. 

The discovery of GRB afterglows at X-ray, optical and radio wavelengths (Costa et 
al 1997, van Paradijs et al 1997, Frail et al 1997) quickly led to a final settlement of the 
issue concerning distance scale: the high redshift associations of GRB-970508 (Metzger et 
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al 1997) and GRB-971214 (Kulkarni et al 1998) indicate that GRB must be of cosmological 
origin, and that its observed fiuence corresponds to a total energy release eo > 10^2-54 gj^gg_ 
In this paper we take Co ~ 4 x 10^^ ergs as representative. Since such a value suggests that 
a process at the end point of stellar evolution (e.g., merging of two neutron stars, collapse 
of a massive star at its evolutionary end point) is responsible for triggering the GRB, it 
is reasonable to assume that the radius of the trigger is Tq ~ 10 km. To the lowest order, 
therefore, one may consider the uniform and simultaneous release of energy to over a radius 
To, and deduce an initial energy density Uo ~ 10^^ ergs cm~^ (hereafter Uo will be expressed 
in such a unit, i.e., f/34). Any internal spatial gradients will after all be rapidly smoothed 
out in the course of the expansion. 

The classic difficulty in explaining a GRB lies not with the availability of a violent 
energy release mechanism, but rather with our ability to account for the basic observed 
properties of the burst. A FB of energy density f/34 = 1 is extremely optically thick to 
Compton- related quantum processes, so that the thermalization timescale is very short. The 
plasma explosion (i.e. 7, e+e~, and baryons) is therefore unobservable electromagnetically, 
and when the plasma expands outwards to become optically thin, its flow turns relativistic 
in situ. The observed non-thermal spectrum and the long duration (1 - 10 s) of a GRB 
implies that most of the energy we detect is from the bulk motion, not the radiation, of the 
outflow (Meszaros and Rees 1993, Paczynski 1993). The conversion back to random energy 
takes place via the extremely rapid particle acceleration at a relativistic shock (Quenby 
and Lieu 1989) when the flow is flnally terminated by the ambient medium, and the energy 
is then radiated away more leisurely. This notion of 'delayed emission' was predicted 
(Meszaros & Rees 1997, Paczynski and Rhoads 1993, Vietri 1997) and is conflrmed by the 
discoveries of afterglows at longer wavelengths. 

The evolution of the FB from creation to observability is, however, easily jeopardized 
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by the presence of trace baryons in the flow, which drastically slow down the bulk motion. 
Indeed, canonical FB models (Meszaros and Rees 1993) presume the mass Mf, of baryons 
in the FB obeys ^ = M^/Mq < 10""^ = C,c (resulting in an outflow that reaches maximum 
bulk Lorentz factor of a few x 10^) although hitherto there does not exist a satisfactory 
explanation of how the requirement may be met. This Letter addresses a new physical 
mechanism which may provide vital clues. 

It has long been recognized (Schwinger 1951) that when a static electric field exceeds 
a critical value E = in^^c' /eh = E^, corresponding to an electron acceleration of ac ~ 2.4 
X 10'^^ cm s~^, it is unstable with respect to pair production. Under the circumstance, a 
virtual e"'"e~ pair is accelerated in opposite directions by the electric field to the speed of 
light within a Compton wavelength. This greatly reduces the probability of annihilation: 
real pair creation may then take place at the energy expense of the field. The process is 
analogous to a lightning discharge, except that it happens in vacuum. The rate of pair 
production per unit volume is given by the formula Schwinger derived: 

d^N a^E^ ^ „2 f-rmrEA _. _i 

} ^ m exp I — I cm ' s (1^ 



dVdt -K^h ^^ "^ \ E 

where a is the fine structure constant and A^ is the number of pairs created. Using this 
formula, we find that for E r^ Ec the 'vacuum breakdown' causes the field to dissipate its 
energy in a timescale of < 10^^^ seconds, resulting in a gamma ray and pair plasma. 

Owing to the strength of E,., this instability has not been realized in the laboratory. 
However, in the context of a GRB, a short-term but intense surface electric field may 
develop as a result of the large radiation pressure of the explosion. Consider a FB of energy 
density U which expands into its immediately vicinity at speed c, the latter having a proton 
number density rii before the disturbance. Assuming that the fraction of U due to radiation 
(i.e. 7, e+e^) is ~ 100 %, the Eddington force on an electron is ~ aJJ {ac being the 
Compton cross section). Although in an extreme Compton limit the prevalence of forward 
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scattering tends to undermine this force, the Eddington value remains correct if the mean 
free path is short enough that repeated Compton events can 're-use' the scattered radiation, 
which is hkely to be the situation just outside a GRB FB. The radiation will accelerate the 
electron away from its neighbouring protons until a counteracting electric field develops 
between them. At this point the separation distance Ar reaches a value / given by: 

aM = ATiUieH (2) 

Non-relativistically the electrons continue to surge forward until Ar is at the maximum 
value of 2/, and subsequently Ar oscillates between zero and maximum at the plasma 



frequency Up = J Aimie^ / rUe while the two layers of charges accelerate outwards. The 
'equilibrium separation' / as depicted in (2) is dependent on n^. Now the baryon density of 
the stationary medium at the surface of FB is unlikely to be orders of magnitude smaller 
than that inside the FB. Assuming that the density ratio is ?] ~ 1 - 10 %, and that the FB 
density within is initially ~ a few x IQ^^^^^ cm^^ (i.e. ^ ~ 0.01 - 0.5 for To = 10 km), we 
shall express rii in units of 10^^ cm^'^, i.e. n^Q. Since the radiation blast comprises photons 
of energy kT ^ meC^ where aT'^ = U, the high energy limit of the Compton cross section 
may be used: 

a, = 8.8 X 10'^' U'jf (Us,) cm^ (3) 

where 

/(f/34) = 1 + 0.04 In t/34 (4) 

/ in (2) is then given explicitly by: 

/ ~ 3 X 10-''n^,'uljiUsi) cm (5) 

Note that because in the case of an initial GRB FB, U is large enough to induce 
'super-Schwinger' acceleration a ^ Qc (i.e. f/34 > 1), the electron will reach speed c 



-6- 



having only traversed a small fraction of /, and relativistic corrections apply to most of the 
motion. Nonetheless there remains the oscillatory nature of Ar, now defined as the charge 
separation measured with respect to our laboratory system S, so pertinent results may still 
be estimated. For instance, at maximum Ar the electrons are at rest with respect to S, and 
the radiation and electric forces are once again given by acU and Airnie^Ar respectively, 
with the latter having exceeded the former. Thus / as given by (5) provides a lower limit to 
the growth of the electric field. The field will oscillate between E = and E > (TcU/e as 
the separation oscillates between Ar = and Ar > /. 

If cXcU/e > Ec the field can, in principle, discharge into pairs once the separation Ar 
is large enough that E ~ Ec. In reality, however, timescale comparisons indicate that the 
field can reach peak strength, with discharge happening efficiently at this point only if 
conditions are favorable (see below). Microscopically the pair production process returns an 
electron to the proton side of the double-layer and creates a positron on the electron side 
(see Figure 1). The pairs are free to escape, as they are created outside the FB, and as as 
we ignore the effects of a magnetic field in this first attempt on the problem. Moreover, the 
discharge is sustained at the FB surface so long as a super-Schwinger surface field continues 
to be regenerated by the radiation pressure of a homogeneously expanding FB. 

To further investigate the Schwinger mechanism, equation (1) may be expressed in 
energy units since the temperature T of the post- discharge e^e~ pairs is obtainable by 
equating the energy density at peak field 

^ ~ ^ r-^y (6) 

with aT"^. The mean energy of each pair, 3kT (^ meC^), is then multiplied by (1) to give a 
peak volume energy loss rate due to discharge: 

^ = 4.7 X 10''ui[f{Uu)]h{U34) ergs cm'^ s-^ (7) 
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where 

9{U3a) = E m-^exp\-7.5l x 10-4f/3/[/(f/34)]-'m (8) 

m=l '^ ^ 

is ~ 1 for f/34 > 10~^ and <^ 1 otherwise (the latter means E < Ec, i.e. the low radiation 
pressure only induces a sub-Schwinger surface field, so no discharge takes place). A simple 
division of (6) by (7) yields the timescale of peak discharge, while the time for the field to 
reach this maximum is ~ l/c. The ratio of these two times is given by: 

^""^"*^"" = 29.0nseUjrh~' (9) 

I formation 

and is > 1 for the parameter regime of concern. The power of pair production at the FB 
surface, de/dt, is the product of -^^ and Airr'^l where r is the FB radius and / is given by 



(5). This yields: 

^^1 21 7 

= 1.8 X lO^'n^y^.U^l [f{Usop9{U,,) ergs s'^ (10) 
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where rio is r in units of 10 km. 



The rate at which electrostatic field energy develops to its peak value at the FB surface 
is given by the product -^ x Anr'^l x y, and may be expressed as: 



de 

'di 



= 5.1 X 10^Vfof/3^[/(f/34)]' ergs s"^ (11) 

field 



where use has been made of (3) and (4). The excess of (11) over (10) is by the same ratio 
as that of the field discharge to formation times given in (9). However, for an efficient 
conversion of the FB energy into surface pairs, one not only must equalize the two rates 
(10) and (11) at peak field, but this rate must also be > the outward energy fiux of the 
explosion itself. For the parameters of the initial FB used here, the total energy fiux across 
a given surface due to the expanding FB is higher than the rates in (10) and (11), and is 
given by: 

= Anr^cU = 3.8 x lO^V^of/34 ergs s"^ (12) 



de 
di 



FB 



In fact, when parity between (11) and (12) is reached, we have the radius independent 
relation: 

which requires an energy density of the surface electric field (6) equal to that of the FB. For 
f/34 ~ 1 the former falls short of the latter by a factor ~ 10. We shall, however, discuss two 
possible scenarios (not exhaustive) under which the rates (10) and (11) may be comparable, 
and be on par with (12) at some early phase of the explosion. The Schwinger mechanism 
may then be important. 

The first scenario maintains our default values of Tq and Uo, but takes advantage of the 
fact that in the foregoing development U refers to the energy density of the FB apparent 
to a laboratory observer at rest w.r.t. the ambient medium, the only exception being (3) 
and (4) where U is the black body energy density at the mean laboratory frequency of the 
radiation blast. As the FB expands a little, and its bulk flow develops a Lorentz factor 
7, volume expansion reduces the co-moving number densities of the FB by 1/7'^, while 
adiabatic cooling scales (Meszaros, Laguna and Rees 1993) the co-moving energy density as 
1/7^ and the energy of a typical photon as I/7. Upon transformation to the laboratory 
system, the length contraction effect, relativistic beaming, and the 'blueshift' of particle 
energies by a factor ~ 7 imply that the photon frequency, number density, and energy 
density are increased from their corresponding co-moving values by a factor of 7, 7^ and 
7® respectively. Thus, at finite 7 the effect of the FB on the ambient medium should be 
calculated using number and energy densities which are higher than their initial values by 
7^ times (e.g. U ~ 7^f^), except ac, which remains unaffected because the photon blueshift 

is counteracted by adiabatic losses. 

_ 1 
The pair efficiency requirement (13) may then be met at 7 > 3[f{U34)]^^U^4^'^. For a 

FB with Uo of f/34 ~ 1 the expansion can easily reach 7 > 3 (at which point > 50 % of f/ is 
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due to radiation). The criterion on ^ for this to happen, ^ < 8.8 x 10^^ = .^o, is consistent 
with an outflow from neutron star merging (Lattimer and Schramm 1976), and is much 
less stringent than that of canonical FB models, since ^/^c ~ 100. On the requirement of 
rapid discharge the ratio in (9), now modifled by the fact that a^ does not scale with 7, 
is ~ 1 for n36 > 1. Under these conditions most of the FB energy e^ will be converted 
to pairs, and will eventually emerge as a GRB. Since the initial FB baryon density is 
> 2.5 X IQ^^ {^ / ^oY'ii cm^^ , and the laboratory density of the FB at the time of discharge is 
7^ times higher, the density ratio across the active front is a reasonable r^ ~ a few %. 

In our second scenario an efficient discharge occurs at the initial explosion 7 = 1. 
A larger value for Uo at the trigger, f/34 > 20, or eo ~ 9 x 10^^ ergs for Tq ~ 10 km, is 
assumed. This will be a limit where (11) is slightly less than (12), so that even though (10) 
can be as large as (11) for nag > 2, only ~ 80 % of to is available to drive a GRB. Under 
optimal conditions, therefore, the burst will have total energy ~ 7 x 10^^ ergs, which may 
be relevant to the more energetic events such as the one detected recently (Kulkarni et 
al 1998) There is no restriction at all on S,'-, but for ^ < 1 the FB is radiation dominated 
and has energy < MqC^. The FB has number density of baryons < 3 x 10^^^ cm~^ and 
kT > 200 MeV, which renders it fortuitously resembling the inner core of a supernova. The 
density ratio rj > 0.7/.^ %. 

Since the region outside the created pairs is essentially the ambient medium which is 
likely to be optically thin to the pairs, a conservative estimate of the baryon contamination 
is obtained by assuming that all the baryons of the original electric double-layer are 
dragged out after discharge. Then the Entropy-per-baryon of the ffow is > f//(nj?72pC^), 
and is > 1/1] times higher than that of the FB itself. The formula provided by Meszaros, 
Laguna and Rees (1993) shows that the ffow will reach a maximum speed given by 
7 > 2.4 X 10^^(eo/4 X lO^'^ergs)^^^'!]^^ . In the ffrst scenario where eo ~ 1 and both rj, ^ are 
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~ 1 %, 7 > a few X 10^, well within the range of values capable of delivering a GRB at the 
terminal relativistic shock. Similarly in the second scenario where eo ~ 20 and r] ^ 1 %, we 
have 7 > 100 if ^ ~ 0.5. 

In conclusion, this Letter demonstrates that a fundamental quantum process which 
takes place at the surface of a powerful celestial explosion has hitherto been ignored. 
The resulting pair outflow can have the right properties to form a GRB at large radii. 
If the proposed mechanism is relevant, further research on the details of its operational 
environment may shed light on the origin of GRBs. 

We thank Martin Rees, Bohdan Paczynski, Stan Woosley, John Doty and Walter Lewin 
for helpful discussions. 
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Figure Caption 

Figure 1. Left: radial outflux of radiation induces a surface electric field (sometimes called 
a Pannekoek-Rosseland field) by accelerating ambient electrons away from protons. Right: 
if this electric field exceeds E^ its discharge will result in the creation of e^e~ pairs outside 
the FB. The pairs can thus expand without encountering more FB matter. 
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